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The chiral geometry of the multiple chiral doublet bands with identical configuration is discussed 
for different triaxial deformation parameters 7 in the particle rotor model with rvhn /2 ® The 

energy spectra, electromagnetic transition probabilities B(M 1) and B(E 2), angular momenta, and 
A-distributions are studied. It is demonstrated that the chirality still remains not only in the yrast 
and yrare bands, but also in the two higher excited bands when 7 deviates from 30°. The chiral 
geometry relies significantly on 7, and the chiral geometry of the two higher excited partner bands 
is not as good as that of the yrast and yrare doublet bands. 
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I. INTRODUCTION 

Spontaneous chiral symmetry breaking phenomenon 
captures significant attention nowadays. Since the chi¬ 
rality of a triaxial atomic nucleus was first predicted in 
1997 by Frauendorf and Meng [1] and identified in 2001 
by Starosta et al. [2], plenty of achievements have been 
made. 

Experimentally, more than 30 candidate chiral nuclei 
have been found at A ~ 80 [3, 4], 100 [5-20], 130 [2, 21- 
46] and 190 [47-51] mass regions. More details see, e.g., 
Ref. [52-57]. 

Theoretically, many approaches have been come up to 
investigate nuclear chirality, such as the particle rotor 
model (PRM) [1, 25, 30, 58-76], the tilted axis crank¬ 
ing model (TAC) [1, 77-79], the tilted axis cranking plus 
random phase approximation (TAC+RPA) [39, 80], the 
interacting boson fermion-fermion model (IBFFM) [37, 
40, 81-84], pair truncated shell model (PTSM) [85-87], 
projected shell model (PSM) [88]. Recently, based on 
the TAC, a microscopical collective Hamiltonian was con¬ 
structed and applied to the unified description of chiral 
vibration and rotation [89]. 

Due to the successful description of nuclear global 
properties and exotic phenomena [90-92], the covariant 
density function theory (CDFT) is introduced to investi¬ 
gate the triaxial deformation of chiral nuclei with various 
configurations. Based on adiabatic and configuration- 
fixed constrain triaxial CDFT, the possible existence of 
multiple chiral doublets (MyD), i.e., more than one pair 
of chiral doublet bands in one single nucleus, was pro¬ 
posed for odd-odd nucleus 106 Rh [93]. Then more MyD 
nuclei were predicted in the other rhodium isotopes [94- 
96] . Later on, the first experimental evidence for the pre¬ 
dicted MyD was reported in 133 Ce [97], and also possibly 
in 107 Ag [98], 

In contrast to the MyD in which different partner 
bands are of distinct triaxial deformations and config¬ 


urations, the MyD phenomenon are expected also with 
identical configuration, i.e., not only the yarst and yrare 
bands but also higher excited bands might be chiral part¬ 
ner bands [99-101]. The energy spectra of MyD are first 
presented with rigid and soft cores [99]. Then the chi¬ 
ral geometry of MyD was examined for configuration 
7r ^u/2®i'/i]] 1 /2 with constant and spin depended variable 
moments of inertia at 7 = 30° [100]. It was confirmed 
that the MyD indeed can exist in identical configuration. 
The chiral geometry of MyD has been further examined 
for the configuration 7r'<79/2 ® at 7 = 30°, 20° [101]. 

Very recently, the evidence of this type of MyD was first 
observed in 103 Rh [102]. This observation shows that 
the chiral geometry in nuclear can be robust against the 
increase of the intrinsic excitation energy. 

As mentioned previous, the MyD with identical con¬ 
figuration has been investigated for the ideal chiral sys¬ 
tem, i.e., one h n / 2 proton particle and one hn/ 2 neutron 
hole coupled to a rigid rotor with triaxial deformation 
7 = 30° [100], and its chiral geometry is examined by an¬ 
alyzing the evolution of angular momentum. Thus, it is 
intriguing and necessary to study the chiral geometry of 
MyD with identical configuration in more general cases 
where 7 deviates 30°. 

In this paper, the physical observables and chiral ge¬ 
ometry of the MyD with identical configuration are dis¬ 
cussed for different triaxial parameters 7 with 7r h\\/ 2 <8> 
in PRM. The paper is organized as follows. In 
Sec. II, the theoretical framework of PRM is briefly in¬ 
troduced. The numerical details are presented in Sec. III. 
In Sec. IV, the obtained energy spectra, electromagnetic 
transition probabilities B(M 1) and B(E 2), angular mo¬ 
menta, as well as AT-distributions are shown and dis¬ 
cussed in details. Finally, the summary is given in Sec. V. 


II. THEORETICAL FRAMEWORK 


- The detailed theoretical framework for PRM can be 

found in Ref. [1, 58, 61, 67, 74]. In this section, for com- 
* Electronic address: qbchen@pku.edu.cn pleteness, some key formulas are presented. The Hamil- 
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where the indices k = 1,2,3 refer to the principal axes 
in body-fixed frame of reference, and 4 , 74, 4 repre¬ 
sent the angular momenta of the total nucleus, core, and 
valence nucleons. The moments of inertia along three 
principal axes depend on triaxial deformation parameter 
7 , Jk = Josin 2 (7 — 2nk/3) [103]. The Hamiltonian of 
intrinsic nucleon H lntr takes the single-j shell Hamilto¬ 
nian [1, 58, 61, 67, 74], 

The eigenstates of PRM are obtained by diagonalizing 
the Hamiltonian (1) at strong coupling basis [58] 
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where | IMK) is the Wigner D function, | k p k n ) is the 
product of the proton and neutron states, is the 

expansion coefficient. The angular momentum projec¬ 
tions onto the quantization axis (3-) in the intrinsic frame 
and the z axis in the laboratory frame are denoted by K 
and M , respectively. The other quantum numbers are 
denoted by a. 

With the obtained eigenstates, the reduced transition 
probability can be calculated according to 
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Here, Q o is intrinsic quadrupole moment; gn, g p , and g n 
are gyromagnetic ratios of rotor, proton, and neutron. 


III. NUMERICAL DETAILS 

In Refs. [1] and [100], the PRM with a particle-like 
/in /2 proton and a hole-like /in / 2 neutron for 7 = 30° 
is respectively used to study the chiral doublet bands 
and the MyD with identical configuration. In Ref. [ 68 ], 
the same system is discussed with different triaxial de¬ 
formation parameters 7 = 30°, 27°, 24°, 21°, 18°, and 
15° to examine the B(M 1) staggering as a fingerprint 
for chiral doublet bands. In this paper, the same model 
will be used to examine the robustness of MyD with 
identical configuration with respect to the triaxiality of 
nucleus. For all the PRM calculations, we employ a 
symmetric particle-hole configuration 7 r/in/ 2 ® Z/ '^'h/ 2 > 
quadrupole deformation with /3 = 0.25, and the mo¬ 
ment of inertia J = 30 MeV _ 1 /i 2 . For the electro¬ 
magnetic transition calculations, the empirical intrinsic 
quadrupole moment Q 0 = (3/\Z5ttRqZ/3) = 3.5 eb, gy¬ 
romagnetic ratios gn = ZfA = 0.44, and g p = 1.21, 
g n = — 0.21 are used [ 68 , 100 ] in accordance with the 
mass region A ~ 130. The numerical details used here 
are in agreement with Ref. [ 68 ] , and only /3 is sightly dif¬ 
ferent (/3 = 0.22 in Ref. [ 68 ]). Thus, the discussions of 
bands 1 and 2 presented in Ref. [ 68 ] still hold true here, 
and we focus on the higher excited bands 3 and 4. 


where a = E or M indicates electric or magnetic transi¬ 
tion, respectively, and A is the rank of electric or magnetic 
transition operator from initial state |i) to final state |f). 
Thus, for reduced electric quadrupole transition proba¬ 
bility 
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IV. RESULTS AND DISCUSSION 
A. Energy spectra 

Fig. 1 shows the energy spectra of four lowest bands 
1, 2, 3 and 4 calculated in PRM when triaxial deforma¬ 
tion parameter 7 varies from 30° to 15°. The discussions 
about the changing of the energy spectra with respect 
to 7 for bands 1 and 2 can be found in Ref. [ 68 ]. For 
7 = 30°, it is difficult to distinguish the doublet bands 
1 and 2 . As 7 deviates from 30°, the energy difference 
between the bands 1 and 2 increases. 

Similar behavior can be also found in doublet bands 
3 and 4. For 7 = 30°, analogous to bands 1 and 2, the 
energy of bands 3 and 4 is identical when 16 < I < 194. 
The energy differences between these two bands vanish, 
indicating that the static chirality appears. In addition, 
the energy differences between neighboring bands at low 
spin region (8 < I < 12 h) are approximately same, which 
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may correspond to the chiral vibration [39, 67]. When 7 
gradually deviates from 30°, the degeneracy at 16 < I < 
19 h is gradually removed. Especially when 7 = 15°, the 
two energy difference between bands 3 and 4 is larger 
than 0.4 MeV, suggesting the absence of static chirality. 
In contrast to the static chirality, the chiral vibration 
remains when 7 decreases. 


—“— Band 1 —■— Band 2 —°—Band 3 —•—Band 4 
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FIG. 1: (Color online) The energy spectra of the four lowest 
bands for the configuration 7 ® with deformation 

parameter /3 = 0.25 and moment of inertia J = 30 MeV _ 1 ft 2 
calculated by PRM with different triaxial deformation param¬ 
eters 7 . Energy of a rigid-rotor has been subtracted from the 
total energies. 


B. Electromagnetic transition probabilities 

Besides the small energy difference, the similarity of 
the electromagnetic transition probabilities and the odd- 
even staggering of B(M 1) are remarkable characteristics 
for chiral doublet bands [36, 37, 60, 64, 68 , 100]. It has 
been demonstrated that for the yrast and yrare bands, 
the B(M 1) staggering is weak in the chiral vibration re¬ 
gion, while it is strong in the static chirality region [ 68 ]. 
Here, we discuss the electromagnetic transition probabil¬ 
ities for bands 3 and 4. The calculated intra- and inter¬ 
band reduced magnetic dipole transitions probabilities 
B(M 1) and electric quadrupole transitions probabilities 
B(E2) of bands 3 and 4 by PRM are respectively shown 
in Fig. 2 and 3 with different triaxial deformation param¬ 
eters 7 . Notice that the intra-band transition probabili¬ 
ties of the doublet bands exhibit similar tendency. The 
same phenomena also take place in inter-band transition 
probabilities. This suggests that the chiral doublet bands 
might remain in bands 3 and 4 when 7 deviates 30°. 

As shown in Fig. 2, for 7 = 30°, when / < 15H, 
the intra-band B(M 1) decreases gradually with spin, 
while the inter-band Ml transitions are forbidden. For 
/ > 15H, a strong staggering between bands 3 and 4 


—•—Band 3—>3 —°—Band 4->4 
—■— Band 3 —>4 —°— Band 4—>3 
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FIG. 2: (Color online) The intra- (band 3 —> band 3 and band 
4 —> band 4) and inter-band (band 3 —> band 4 and band 4 
—» band 3) reduced magnetic dipole transitions probabilities 
B(M1) of bands 3 and 4 calculated by PRM with different 
triaxial deformation parameters 7 . 


can be clearly seen, which means there is an alternative 
suppression between the intra- and inter- transition as 
spin increases. This has been demonstrated in Ref. [100]. 
However, as 7 deviates from 30°, two obvious changes 
happen. One is that the staggering becomes weaker and 
weaker at high spin region when 7 varies from 27° to 15°; 
the other is that the suppression of inter-band transitions 
is gradually relieved at low spin region. Compared with 
the first partner bands 1 and 2 , for which the staggering 
when 7 = 15° completely disappears [ 68 ], the staggering 
of the second partner bands 3 and 4 still exists at the 
same spin internal. 

Fig. 3 shows the electric quadrupole transitions prob¬ 
abilities B(E 2) between the second partner bands. It 
can be seen that for 7 = 30°, when I < 15/i, the intra¬ 
band E2 transitions are forbidden, and the inter-band 
ones are allowed. When I > 15/i, on the contrary, the 
intra-band E2 transition are allowed, and the inter-band 
E2 transitions are forbidden. As 7 deviates from 30°, 
the probabilities of both intra- and inter-transitions for 
I > 20 h are similar, which indicates that B{E2) is irrel¬ 
evant to 7 at high spin. At low spin (/ < 20 A), the sit¬ 
uation is inverse: the suppressed inter-bands transitions 
are gradually allowed, yet intra-bands are gradually for¬ 
bidden. The similar conclusions for bands 1 and 2 has 
been derived in Ref. [ 68 ]. 

Therefore, from the discussions above, it can be 
found that for / > 18 h, the odd-even staggering of 
B(M 1) will lead to the staggering of B(M1)/B(E2) and 
B(Ml)i n /B(Ml) out . When 7 gradually decreases, the 
amplitude of B(M1)/B(E2) and B(Ml)i n /B(Ml) out 
staggering will decline because of the weakening of 
B(M 1) staggering. 
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FIG. 3: (Color online) The intra- (band 3 —> band 3 and band 
4 —> band 4) and inter-band (band 3 —» band 4 and band 4 —>• 
band 3) reduced electric quadrupole transitions probabilities 
B(E2) of bands 3 and 4 calculated by PRM with different 
triaxial deformation parameters 7 . 


C. Angular momenta 


An insight into chiral geometry will be revealed when 
it comes to the rms values of the angular momenta com¬ 
ponents of the core R k = (R^) 1 / 2 , the valence proton 
Jpk = (jp fc ) 1/2 5 and the valence neutron J nk = (j 2 fe ) 1/2 
(k = 1,2,3). The results for bands 1 and 2 have already 
shown in Ref. [ 68 ] , where the chiral geometry for the dou¬ 
blet bands is discussed. The obtained results for bands 3 
and 4 at triaxial deformation parameter 7 = 30°, 21°, 15° 
are shown in Fig. 4, in which l,i,s correspond to long, 
intermediate and short axes. It is shown that for both 
partner bands, the core angular momentum mainly aligns 
along the Taxis due to its largest moment of inertia. The 
contribution of the valence proton particle mainly along 
the s-axis, and that of the valence neutron hole mainly 
aligns along the Taxis, for such orientations are favored 
by their interaction with the triaxial core [1]. With the 
total angular momentum increasing, R increases gradu¬ 
ally, J n and J p move gradually toward the Taxis due to 
the Coriolis interaction, and the three angular momenta 
together form the chiral geometry of aplanar rotation. 

From the left panel, it can be easily seen that when 7 = 
30°, the angular momenta of bands 3 and 4 are almost 
same, except at band head and / = 20 H. They have 
the similar orientation. However, as 7 gradually deviates 
from 30°, this similarity disappears gradually. When 7 = 
21 °, the differences between the partner bands become 
larger than those for 7 = 30°, yet at I = 20,21,23 H, 
two bands are still similar. As 7 continually decreases, 
such as 7 = 15°, the orientation of angular momenta 
at 20 < I < 24 H of the core, valence proton and neutron 
are different, corresponding to the large energy difference 


between bands 3 and 4 shown in Fig. 1. 
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FIG. 4: (Color online) The expectation values of angular mo¬ 
menta of the core (R k = \J(R\)), the valence proton par¬ 
ticle (Jpk = \J(Jpk)), an d the valence neutron hole (J n t = 

\j(Jnk)) along the intermediate (i), short(s), and long (Z) axes 
with the triaxial deformation parameter 7 = 30°, 21°, 15°. 


D. iv-distributions 

To better understand the mechanism of chiral geom¬ 
etry, the IT-distributions, which reveal the probability 
distributions of the projection on three principle axes of 
the total angular momenta, are studied. 

As discussed above, at low spin, the valence proton 
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FIG. 5: (Color online) Diagrammatic sketch of chiral vibra¬ 
tion of total angular momentum. 


particle and neutron hole contribute most to the total an¬ 
gular momentum. Thus the angular momentum mainly 
lies in s-l plane, as diagrammatically shown in Fig. 5. 
Due to the effect of quantal fluctuation, the total angular 
momentum will vibrate through s-l plane. In this case, 
the chiral vibration will appear since the total angular 
momentum oscillates between the left-handed and right- 
handed systems [39, 67]. To illustrate the vibrational 
character of this dynamical chirality procedure, Fig. 6 
shows the A-distribution along the i-axis of four lowest 
energy bands for 7 = 30° at I = 10ft. The A-distribution 
of bands 1 and 2 have been shown in Ref. [ 68 ] , where it 
has been pointed out that bands 1 and 2 correspond to 
the zero- and one- phonon states with symmetric and an¬ 
tisymmetric wave functions on the i-axis. In the same 
way, bands 3 and 4 represent two- and three-phonon 
states with symmetric and antisymmetric wave functions. 
The chiral vibration character exhibits at low spin region 
can explain the similar spacing between the neighboring 
bands at the beginning of the MyD as shown in Fig. 1. 

The obtained results for bands 3 and 4 at triaxial defor¬ 
mation parameter 7 = 30°,21°, 15° are shown in Fig. 7. 
In the following discussions, we use A s , Ki , and A/ to 
represent the distributions on the short, intermediate, 
and long axis. 

For 7 = 30°, at higher spin region, such as I = 18, 22ft, 
the A-distributions of the partner bands turn to be simi¬ 
lar, implying the alteration from chiral vibration to static 
chirality [ 100 ]. 

As 7 decreases, e.g., 7 = 21°, the A-distribution of 
bands 3 and 4 when I < 20ft is similar to that when 7 = 
30°. But another mode of new type of chiral vibration is 
found when / > 20ft. The maximum K s distribution at 
that internal appears near K s = 0 for band 3 and band 
4, which corresponds to a chiral vibration through the l-i 
plane. 

When 7 = 15°, two modes of chiral vibration remain at 
low and high spin, but the pure static chirality is hardly 
distinguishable. Thus, for spin 16 < I < 20ft, chiral 
vibration is mixed with static chirality, which may char¬ 
acterize energy spectra by the removal of degeneracy at 
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FIG. 6 : (Color online) The A-distribution of four lowest 
bands 1, 2, 3, and 4 for 7 = 30° when I = 10ft. 


16 < I < 20 ft. 


E. Energy spectra for higher excited states 

The discussion above gives the results of the physical 
observables and chiral geometry of the MyD with identi¬ 
cal configuration for different triaxial parameters 7 with 
7 rftii /2 ® vh^ i/ 2 . It has been shown that the chirality 
still remains not only in the yrast and yrare bands, but 
also in the two higher excited bands when 7 deviates 
from 30°. Further consideration of higher excited bands 
seems interesting. In Ref. [89, 100], the energy spec¬ 
tra of higher excited bands 5 and 6 have been shown at 
7 = 30°. It is found that bands 5 and 6 might also be a 
pair of chiral doublet bands. In Fig. 8 , the energy states 
as well as their excitation energies with respect to the 
yrast states obtained by PRM with triaxial deformation 
parameters 7 = 30° are shown as functions of spin. For 
each spin, the lowest twenty states are shown. It can be 
seen that in high excited bands, whose excitation energies 
are larger than 2.0 MeV, the situation turns to be com¬ 
plex. The level density is rather high and neighboring 
bands are too close to be distinguished. Thus it would 
be very difficult, if not impossible, to select the chiral 
doublet bands [100]. In fact, these states would fall into 
the regime of quasi-continuum spectra and the strong in¬ 
teractions between them evoke the emergence of rotation 
damping [104, 105]. This complexity can be explained by 
the PRM Hamiltonian (1), where the single particles, the 
rotor, and their coupling Hamiltonian are involved. If the 
single particles were not considered, it would turn to be 
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FIG. 7: (Color online) The probabilities distributions of the 
projection of total angular momenta of bands 3 and 4 along 
the intermediate, short, and long axes when the parameter 
7 = 30°,21°,15°. 


a triaxial wobbling rotor Hamiltonian [106]. At the spin 
region (8 ~ 20h) discussed above, the wobbling rotor has 


low angular momentum 10 A), which corresponds to a 
complicated wobbling energy spectrum [106-109]. With 
the coupling of the single particles, the occupation of pro¬ 
ton and neutron at different orbitals further complicates 
the wobbling energy spectra of triaxial rotor. 

In addition, the energy spectra of higher excited states 
with triaxiality parameters smaller than 7 = 30° are also 
examined. We note that the triaxial deformation has 
larger influence on the low lying states than high ones. 
When the triaxial deformation parameter varies from 30° 
to 15°, the complexity of the high energy states still re¬ 
main. Anyway, more efforts need to be paid to investigate 
the higher excited states. 


V. SUMMARY 

In summary, the chiral geometry of the M%D with 
identical configuration is discussed for different triaxial 
deformation parameters 7 with nhn /2 0 in par¬ 

ticle rotor model. The energy spectra, electromagnetic 
transition probabilities B(M 1) and B(E 2), angular mo¬ 
menta, and ^-distributions are studied. The calculated 
results strongly suggest the chirality still remains in 
not only the yarst and yrare bands but also two higher 
excited bands when 7 deviates from 30°. The chiral 
geometry relies on 7 , for the chiral vibration and static 
chirality exist in different spin internal when 7 changes. 
In addition, the chiral geometry of the higher energy 
partner bands is not as obvious as the ones with the 
lowest energy. 
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